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We present two schemes for driving Raman transitions between the ground state hyperfine man- 
ifolds of a single atom trapped within a high-finesse optical cavity. In both schemes, the Raman 
coupling is generated by standing-wave fields inside the cavity, thus circumventing the optical access 
limitations that free-space Raman schemes must face in a cavity system. These cavity-based Raman 
schemes can be used to coherently manipulate both the internal and the motional degrees of free- 
dom of the atom, and thus provide powerful tools for studying cavity quantum electrodynamics. We 
give a detailed theoretical analysis of each scheme, both for a three-level atom and for a multi-level 
00 \ cesium atom. In addition, we show how these Raman schemes can be used to cool the axial motion 

c 2 ^ , of the atom to the quantum ground state, and we perform computer simulations of the cooling 

^^ ' process. 

PACS numbers: 32.80.Qk, 32.80.Lg, 32.80.Pj 

\0 \ I. INTRODUCTION 

r^ ■ Systems consisting of a single atom coupled to a high-finesse optical cavity are of fundamental importance to 

Q^' quantum optics and quantum information science. Such cavity QED systems have been experimentally implemented 
using neutral atoms U ] , B , S 0, S Q and ions [3, W\, and have been the subject of numerous theoretical studies 
g.[lQ.O . [l2 .[lall4l.lla.[l6l|. In particular, such s yst ems play a key role in proposals for scalable quantum computation 
ITI . [l8j and distributed quantum networks [l^ [20| . An important requirement for many of these proposals is the 
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ability to coherently control the internal and motional degrees of freedom of the trapped atom, and Raman transitions 
provide the means for meeting this requirement. 

Raman transitions are powerful tools that have diverse applications in atomic physics, including spectroscopy [2l| . 
'T'* ] precision measurement [22, [23[ , and coherent state manipulation [2J] , and have been used to coherently control the 
J_^ ,. motional degrees of freedom of trapped ions [2^ [2y| and of neutral atoms trapped in optical lattices [23, [2^, l2g |. 
_^ Until recently, however, Raman transitions had not been incorporated into cavity QED. The practical challenge of 
^\^ implementing a free-space Raman scheme in cavity QED lies in the presence of the cavity itself, which offers limited 
^^ optical access to the atom within, especially in the strong coupling regime [0] . As first proposed in [S^ , these optical 

I [ \ access limitations can be circumvented by implementing a cavity-based Raman scheme, in which the Raman coupling 

f— ^ ■ is generated by standing-wave fields inside the cavity. 

OQ ' In this paper, we provide the theoretical background behind two such cavity-based schemes for driving Raman 

f^ . transitions between the hyperfine ground state manifolds of a single atom optically trapped within a high-finesse 
cavity. Both schemes have been recently implemented and validated experimentally, and have been used to extend 
the trapping lifetime of a near-resonantly driven atom [SO] , to cool the axial motion of an atom to the quantum ground 
state of its trapping potential [3r| , and to optically pump an atom into a specific Zeeman state [32] . 
J^ I The paper is organized as follows. In section|TTl we describe how a single atom is trapped inside the cavity by means 

of an optical dipole trap. In section IIIIl we present the two schemes for driving Raman transitions in the trapped 
atom. We treat the atom using a three-level model, and show that the internal and motional degrees of freedom of 
the atom can be described using an effective Hamiltonian that has the same form for both schemes. In section HVl we 
quantize the axial motion of the atom, and show how the Raman couplings allow one to drive transitions that change 
the vibrational quantum number for motion along the cavity axis. We present both analytic results, which apply to 
atoms that are sufficiently cold, and numerical results, which apply to atoms of arbitrary temperature. In section FVl 
we describe how the Raman couplings are modified when we take into account the multiplicity of levels in a physically 
realistic cesium atom. We show that the Raman schemes drive transitions between individual Zeeman states in the 
two hyperfine ground state manifolds, and we derive the Rabi frequencies for these Zeeman transitions. Finally, in 
section rvTl we show how Raman transitions can be used to cool the axial motion of the atom to the quantum ground 
state, and we present computer simulations of the cooling process. 



II. TRAPPING AN ATOM INSIDE THE CAVITY 

The optical cavity we will be considering consists of two symmetric mirrors separated by a distance L (details of the 
cavity are given in Appendix]^. An atom is trapped inside the cavity by means of a far off-resonance trap (FORT), 
which is created by driving one of the cavity modes with light that is red-detuned from a dipole transition in the 
atom. The red-detuned light forms a standing wave inside the cavity, and the coupling of the atom to the light causes 
it to be attracted to the points of maximum intensity in the standing wave. 

The operation of the FORT can be understood by considering a simple two-level model in which the atom has a 
single ground state g and a single excited state e. Let We denote the frequency of the g ~ e transition, and let ujp 
denote the frequency of the FORT light. We will assume that the FORT resonantly drives mode np of the cavity, so 
wf = 'i.'KVFSRnF^ where vfsr = 1/2L is the free spectral range. The Hamiltonian for the system is 

H = uje\e){e\ + {Vt + 9)) cos uj Ft, (1) 

where 

Cl = 9F^F{r^\g){e\. (2) 

Here f is the position of the atom, V'-F(r) is a dimensionless quantity that characterizes the shape of the cavity mode 
(see Appendix [X| , and VIf is the Rabi frequency of the light at a point of maximum intensity. We can simplify this 
Hamiltonian by making the rotating wave approximation and then performing a unitary transformation to eliminate 
the time dependence: 

H^-^F\e){e\ + \{h + ^^), (3) 

where Aj? = ojl — ijJe is the detuning of the FORT from the atom (note that because the FORT is red-detuned, 
/S.F < 0). In the limit that the FORT is far-detuned (|Af| ^ ^_f), we can adiabatically eliminate the excited state 
(see [3J|) to obtain an effective Hamiltonian for the ground state: 

HE = ^^(l(l^^U{r)\g){gl (4) 

where 

C/(r) = -C/F|V'i=^(r)|2 (5) 



describes a trapping potential with depth Uf = il|,/4|Ai?|- For the experiments described in [30|, |3l|, |33] , the power 
in the FORT beam is set such that Uf — (27r)(50 MHz). To describe the shape of the potential, it is convenient to use 
a cylindrical coordinate system centered on the cavity axis: we will let z and p denote the axial and radial coordinates 
of the atom, where the cavity mirrors are located at z = and z = L. Using equation (jAip to substitute for the mode 
shape ipF{r) in this coordinate system, we find that 

U{r) = -Uf e-2p'/»F gin^ kpz, (6) 

where kp = Trnp/L is the wavenumber for the FORT mode np- The minima of the potential occur at the points 
p — 0, z — Zr, where Zr is defined such that 

kFZr^Tr{r + l/2). (7) 

Since < Zr < L, we find that r = 1, ...,np — 1; thus, there are np distinct FORT wells in which an atom can be 
trapped. Let us assume that the atom is trapped in FORT well r. It is convenient to define a coordinate x = z — Zr 
that gives the axial displacement of the atom from the potential minimum of this well. We can then express the 
trapping potential as 

U{r) = -Up e-^p'/"-! cos^ kpx. (8) 

Near the bottom of the well the potential can be approximated as harmonic, 

U{r) ~ -Up + -mulp' + ^rnivlx^, (9) 



where ujr and Wa, the radial and axial vibrational frequencies, are given by 

-muol = 2Uf/w%, -mujl = Upkl- (10) 

The corresponding periods 2Tr/uJr and l-n jiOa characterize the timescales for radial and axial motion. For the exper- 
iments described in [Sfl ISlI, l32|, the vibrational frequencies are uor — (27r)(5kHz) and Wa — (27r)(500kHz), so the 
timescale for radial motion is much longer than the timescale for axial motion. We will be interested in describing the 
evolution of system over timescales that are short compared to the timescale for radial motion, but not necessarily 
short compared to the timescale for axial motion. For such timescales we can view the atom as being radially sta- 
tionary, and take p to be a constant parameter that enters into the potential for the axial motion. Thus, dropping a 
constant term, we can express the potential for the axial motion as 

C/(f) = t/pSin^/c^a;, (11) 

where Up = Up e"^'' ^^^ is the axial trap depth at radial coordinate p. 

III. RAMAN COUPLING FOR A THREE-LEVEL MODEL 

A. Effective Hamiltonian 

To show how a Raman coupling can be generated in a trapped atom, let us now consider a three-level model 
in which the atom has two ground states a and 6, which correspond to the ground state hyperfine manifolds of a 
multi- level atom, and a single excited state e. The excited state has energy We, and the ground states a and b have 
energies —Ahf/^, and +Ahf/'2, where Ahf is the ground state hyperfine splitting. The Hamiltonian for the atom is 

Ho=u;,\e){e\ + -AHFa,, (12) 

where 

a,^\b){b\-\a){a\. (13) 

We can generate a Raman coupling between the two ground states by driving one of the cavity modes with a 
pair of beams that are tuned into Raman resonance with the ground state hyperfine splitting of the atom. Let us 
denote the optical frequencies of these beams by a;± — lol ± Sj^/2^ where lo^ is the average frequency and Sr is the 
frequency difference (see Figure [T]). Also, let us define a parameter S = Sn — Ajjf that gives the Raman detuning of 
the beams; we will assume that the beams are tuned close to Raman resonance, so \S\ <ti Ahf- The beams generate 
standing-wave fields inside the cavity, and the coupling of the atom to these fields is described by the Hamiltonian 

Hr = {Cl+ + ri\_)cosu;+t+ {ft- +n^_)cosuj_t, (14) 

where 

n± = fi± v(^) A (15) 



and 



A^i\a) + \b)){e\ (16) 



is an atomic lowering operator. Here r is the position of the atom, tp{r) is a dimensionless quantity that characterizes 
the shape of the driven mode (see Appendix [X)l . and fl± arc the Rabi frequencies of the fields at a point of maximum 
intensity. For simplicity, we have assumed that the a — e and b — e transitions couple to the light fields with equal 
strength. 

The total Hamiltonian for the system is H — Hq + Hr. We can simplify this Hamiltonian by making the rotating 
wave approximation and then performing a unitary transformation: 

H = -A\e){e\+'^AHF^z+B + B\ (17) 



where 



B = ]-{h+ e*'^"*/^ + r2_ e-''*«*/2)^ (^g) 



and /S. = lol — Lde describes the overah detuning of the optical fields from the excited state. We will assume that 
the fields are far-detuned from the atom (|A| 3> ^±), so we can further simplify the Hamiltonian by adiabatically 
eliminating the excited state to obtain an effective Hamiltonian for the ground states: 

He = ]^^HFCT, + ^BB\ (19) 

= \^HF(J, - Ve |^(r)|2 - Ve 1^(7=01' o^ + ^E |^(r)|2 cosSrI + Qe |^(r)|2 a^ cosSnt, (20) 



where 



and 



VE = inl + nl)/A\A\, nE = ^+^-/2A, (21) 



a^ = \a){b\ + \b){a\. (22) 

We will assume that the optical fields are weak enough that Ve, ^e ^ ^hf, so the first term of He dominates and 
to a good approximation the eigenstates of He are |a) and |6). The second term of He describes a state- independent 
level shift, which is analogous to the FORT potential we derived in the previous section. The third term of He gives 
a state-dependent correction to the level shift, which is of order Ve/Ahf and may therefore be neglected. The fourth 
term of He describes a modulation of the state- independent level shift at frequency Sr. Because we have assumed 
that the system is tuned near to Raman resonance, Sr is of the same order as the hyperfine splitting Ahf- For 
cesium, the atom used in the experiments described in [30, ISlli 123 1 ^^^ hyperfine splitting is Ahf = (27r)(9.2 GHz), 
which is much larger than the harmonic frequencies LUr and tUa that characterize the timescales for atomic motion. 
Thus, over the motional timescales the fourth term of He averages to zero and may also be neglected. After making 
these approximations, we are left with 

He = \Ahfo, - Ve \^{^? + ^e \^{^fo, cosSrI. (23) 

We can further simplify the effective Hamiltonian by making the rotating wave approximation and then performing 
a unitary transformation to eliminate the time-dependence: 

He = -f ^- - Ve IVX^H' + 1^e |^(r)P a^. (24) 

This Hamiltonian describes an effective two-level atom with ground state a and excited state b, which is driven by a 
classical field with Rabi frequency flE |'/'('^)P and detuning 6. We will now use this effective Hamiltonian to describe 
two schemes for driving Raman transitions in a trapped atom. 

B. FORT-Raman configuration 

In the first scheme, which we will call the FORT-Raman configuration, the FORT itself forms one leg of a Raman 
pair. To form the other leg of the pair we add a much weaker beam, which we will call the Raman beam, that drives 
the same cavity mode as the FORT but is detuned from the cavity resonance by 6r. This configuration was first 
proposed in [SJ] , and formed the basis of the optical pumping scheme described in [33] and of the cooling scheme 
used in [30| to extend the lifetime of a trapped atom. Let us denote the optical frequencies of the FORT and Raman 
beams by wp and lur ~ up + 5r, and their maximum Rabi frequencies inside the cavity by O^ and ^r. We can then 
apply the results of section HH Al to obtain an effective Hamiltonian He that describes the FORT-Raman pair. We 
will assume that the Raman beam is blue-detuned from the FORT {5r > 0), so we take the J7+ field to be the Raman 
field and the fi- field to be the FORT field; if the Raman beam is red-detuned from the FORT {5r < 0) then we 
reverse these identifications. 

As was discussed in section [Hi we treat the atom as being radially stationary, and consider only the axial motion. 
Thus, the total Hamiltonian for the atom is given by 

H^^+He. (25) 
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FIG. 1: Level diagram for the three-level atom. Shown are the ground states a and b, the excited state e, and the pair of 
optical fields at frequencies a;+ and lo^ . 

Here p is the momentum of the atom in the axial direction, and He is given by equation ([24]) with i^E = ^f^r/'^^f 
and Ve = Uf + Ur, where Up = ^'f/^\^f\ and Ufi = r^^/4|Ai;'| are the maximum level shifts due to the FORT and 
Raman fields individually, and A^ = LOp — ^e is the overall detuning of the FORT and Raman beams from the atom. 
It is convenient to express H as H = H^xt + Hint, where 



He.t^-^-VE\i^F{r^? 



describes the axial motion of the atom, and 



A 1 

Hint = ^-^^z + -^ ^B \i'F{r)^ Ox 



(26) 



(27) 



describes the internal state of the atom. Typically the powers of the FORT and Raman beams are such that ^e — 
(27r)(200kHz) and Up ^ (27r)(50MHz), and for these values Uh/Up ^ [np/Upf ~ 2 x 10"^ Thus, we can neglect 
the level shift of the Raman beams and approximate Hext as 



Zm 



(28) 



where U{f) — —Up\i{jp{r)\'^ is the FORT trapping potential. As in section [TTl we will assume that the atom is 
trapped in FORT well r, and define a coordinate x ~ z — Zr that gives the axial displacement of the atom from the 
well minimum. Note that 



\^p{r)\^ = e- 



"" sin^ kpz = e' 



"" cos^ kpx, 



(29) 



where we have substituted for the FORT mode shape tjjp{r) using equation (jAip and for kpz^ using equation ([7]). 
Thus, we can express the internal Hamiltonian as 



Hint 



2^z + -;;^pCOS^kpxax, 



(30) 



where Vtp = Q,e e "^p l^f is the effective Rabi frequency for an atom at radial coordinate p. 



C. Raman- Raman configuration 



For the second configuration, which we will call the Raman-Raman configuration, the FORT drives mode np of the 
cavity at frequency cji?, and a pair of Raman beams drives mode n/j of the cavity at frequency lofi. This configuration 
was used to perform the ground-state cooling described in [31|. We will assume that the Raman beams have equal 
powers and are tuned symmetrically about the cavity resonance, so the frequencies of these beams can be expressed 
as wpi ± 5iil2. Let Q,p denote the maximum Rabi frequency for the FORT beam, and let fi/j denote the maximum 



Rabi frequency for one of the Raman beams. Using the results of section UlI Al we can describe the couphng of the 
atom to the pair of Raman fields in terms of an effective Hamiltonian He- 

As was discussed in section |lll we treat the atom as being radially stationary, and consider only the axial motion. 
Thus, the total Hamiltonian for the system is 

H=^ + U{r) + HE. (31) 

The first term describes the kinetic energy of the atom due to axial motion, and the second term describes the FORT 
potential U{r) = -Up \'il)F{r)?, where Up = rj|/4|AF| is the FORT depth. The third term is given by equation ^ 
with Ve ~ 2[/fl, Ur — rj|j,/4|Afl|, and Q,e — ri^/2Afl, where A/j = u)b. — u)^ is the overall detuning of the Raman 
pair from the atom. It is convenient to express H as H = H^xt + Hint , where 

He.t = ^ + U{r) - VE\^R{r)\^ (32) 

describes the motion of the atom, and 

Hint = --x<Jz + -x^E\i'R{r)\^ ly^ (33) 

describes the internal state of the atom. 

Note that because the FORT and Raman beams drive different cavity modes, the registration of the standing waves 
corresponding to the two beams depends on the axial position of the atom, and therefore on the particular FORT 
well in which the atom is trapped. This is in contrast to the FORT-Raman configuration, for which the FORT and 
Raman standing waves are always perfectly registered. One consequence of the well-dependence of the registration is 
that the level shift Ve due to the Raman pair distorts different FORT wells in different ways. We calculate this effect 
in Appendix iBl but for now we note that for the typical parameters VLe — (27r)(200kHz), Up — (27r)(50MHz) the 
ratio of the level shifts due to the FORT and Raman beams is Ur/Uf = ^e/Uf — 4 x 10~^. Thus Ve <C Uf, so we 
can neglect the Ve term and approximate Hf^xt as 

H,xt ^^ + U{r). (34) 

Zra 

As in section |TT1 we will assume that the atom is trapped in FORT well r, and define a coordinate x = z — Zr that 
gives the axial displacement of the atom from the well minimum. Note that 

\i^R{r)? = e^^"'/""" sin^ Urz = e^^pV^'^j cos^{kRX + a), (35) 

where a — {Ur — kF)zr is the phase difference between the FORT and Raman beams at the bottom of FORT well 
r. Here we have substituted for the FORT mode shape tpF{r) using equation (jAip and for kpZr using equation ([7]). 
We will assume that the FORT and Raman beams drive nearby modes of the cavity, so |nfl — npl ^ riR, np- In this 
limit {kR — kp)x <SC 1 and wr^ ~ wp, so we can approximate equation (j35[) as 

|V;fl(f)|2 = e^^"'/'"^ cos^ikpx + a). (36) 

Thus, we can express the internal Hamiltonian as 

A 1 

Hint = - xo-^ + X ^p cos^ {kpx + a)ax, (37) 

n 2 / 2 

where flp = Up e~ ' ""^ is the effective Rabi frequency for an atom at radial coordinate p. 

D. Summary of Raman schemes 

We have shown that that for both the FORT-Raman and the Raman-Raman configurations the Hamiltonian has 
the form H = H^xt + ^mt j where 

Hext = 7, ^Up sin^ kpx, (38) 

Zm 

Hint = --az + -i^p cos^ {kpx + a) ax- (39) 



For the FORT-Raman configuration a = 0, and for the Raman-Raman configuration a = (fc/j — kF)zr for an atom 
trapped in FORT well r. The Hamiltonian Hf,j.t describes the motion of the atom in the axial potential and is 
independent of the atomic state, and the Hamiltonian Hi„t describes a Raman coupling between the ground states 
that depends on the axial position of the atom. 

For either the FORT-Raman or the Raman-Raman configuration, the Raman coupling described by Hint can be 
turned off by turning off one of the beams in the Raman pair (note that for FORT-Raman configuration the FORT 
beam must always be on in order to maintain the trapping potential, so for this configuration one must turn off the 
Raman beam). Alternatively, one can keep the beams in the Raman pair on at all times and turn off the Raman 
coupling by tuning the pair out of Raman resonance {\6\ ^ fip). As we have seen, there is a small level shift due to 
the Raman beams, which we neglected when writing down the above expression for Hg^ti and this detuning-based 
method has the advantage that these level shifts are always present regardless of whether the Raman coupling is on 
or off. With the first method, these level shifts cause a slight change in the trapping potential whenever the Raman 
coupling is turned on or off, which could potentially heat the atom or cause other problems. 

IV. QUANTIZATION OF AXIAL MOTION 

For many applications, such as Raman sideband cooling, it is necessary to quantize the axial motion; that is, to 
treat the axial position x and momentum p as quantum operators. We first show how this is achieved for cold atoms, 
and then discuss some numerical results that apply to atoms of arbitrary temperature. 

A. Approximate form of the Hamiltonian for cold atoms 

For cold atoms the axial trapping potential is nearly harmonic, where the harmonic frequency uj is given by 

w = (2C/p/m)i/2 kp = uja e-"'/"'^. (40) 

We can quantize the axial motion by introducing phonon creation and annihilation operators operators h^ and 5, 
which are related to x and p by 

x= (2TOw)-i/2(6-H6t), p^-i{muj/2)^'^{b~b^). (41) 

From these relations, we find that 

kFX^T^{b + h^), (42) 

where 77, the Lamb-Dicke parameter, is given by 77 = {2mijj)^^i'^ujF- Note that because uj depends on the radial 
coordinate p, the Lamb-Dicke parameter also depends on p. 

If the atoms are sufficiently cold, we can obtain a reasonable approximation to H — Hint + Hext by expanding in 
ry and retaining terms only up to second order; from equations (|38p . ([M)) . and (|42p. we find that 

H,xt = Loil/2 + b^b)-Loir,yU)ib + b^)\ (43) 

H,„t = --a^ + -np{il/2){l + cos2a)-r]{b + b^)sin2a-ri'^{b + b'^fcoa2a)ax- (44) 

It is convenient to form a basis of states {|a, n), \b,n)} for the system by taking tensor products of the internal states 
I a) and |6) with the motional Fock states {|n)}. To order jf these product states are eigenstates of H^xt, where pairs 
of states I a, 71), \b,n) with the same vibrational quantum number n are degenerate and have energy 

En = {a, n\Hext\a, n) = (6, n\H,xt\h, n) = lu{\/2 + n) ~ Lo{rf/A){\ + 2n + 2^2). (45) 

The Raman coupling described by Hint drives transitions between the product states. By taking matrix elements 
of the Raman coupling, we find that state |a,7i) is coupled to states |6, ti), |fe, 7i± 1), and |&, 7i± 2), where the Rabi 
frequencies for these transitions are given by 

^n^nl^p = 1/2 +(1/2 -772(272+1)) COS 2a, (46) 

^n^ri±i/^p = -7]Vn± 1 sin2Q;, (47) 

^n^n±2/^p = -if^n ± l^/n ± 2 COS 2a. (48) 
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Note that An = ±1 transitions are suppressed relative to An = transitions by ^ ily/n, and An = ±2 transitions 
are suppressed relative to An = transitions by ~ rj^ n. To resonantly drive the n ^ n transition we set 5 = 0, and 
to resonantly drive the n — > n ± 1 and n — > n ± 2 transitions we set 6 ~ 5n^n±\ and 8 = Sn^n±2, where 



+n±l 



En±l — En — ±(5„, Sn-tn±2 = En±2 — En — ±2J„, (49) 



and Sn = to ~ ifuj n. For a harmonic trap Sn = tu and the frequencies of these transitions arc independent of n, but 
because the FORT is shallower than its harmonic approximation, Sn decreases with increasing n. 

Recall that for the FORT-Raman configuration a = 0, whereas for the Raman-Raman configuration the value 
of a depends on the FORT well in which the atom is trapped. Thus, in the FORT-Raman configuration the Rabi 
frequencies are the same for all the FORT wells, whereas in the Raman-Raman configuration the Rabi frequencies vary 
from well to well. Also, note that in the FORT-Raman configuration the An = ±1 transitions are always forbidden. 
This follows from symmetry considerations: since the trapping potential is symmetric under x — > —x the motional 
eigenstates are also parity eigenstates, and since the Raman coupling is symmetric under a; — ^ — a; it cannot couple 
an even parity state to an odd parity state. 

B. Numerical results 

The harmonic approximation we described in the previous section only applies to Fock states \n) for which rj\/n <^ 1. 
For higher- lying energy levels, we can calculate the Rabi frequencies and detunings for the various motional transitions 
by numerically solving the time-independent Schrodinger equation for Hf^^t- This provides us with a set of motional 
eigenstates {|V'n)} and eigenvalues {En}- Using the motional eigenstates, we can take matrix elements of the Raman 
coupling described by Hint to calculate the Rabi frequencies for different motional transitions: 

VLn^n/^p = {lpn\cOS^{kFX + a)\'ll)n), (50) 

r2„^„±i/rip = {ipn±i\cos^{kFX + a)\i)n), (51) 

^n^n±2l^p = {lpn±2\cOS^{kFX + a)\i!n)- (52) 

Note that Hint can also drive higher-order n-changing transitions, but the matrix elements for these transitions are 
quite small and we will not consider them here. From the energy eigenvalues, we can determine the detunings for the 
n ^ n ± 1 and n — *■ n ± 2 transitions: 

5n~^n±l — En±l — En, <5n^n±2 — En±2 — En- (53) 

The numerically-determined Rabi frequencies and detunings are shown in Figure [2] for a = and a = 7r/2, and in 
Figure [3] for a = 7r/4. Note that for a = and a = 7r/2 the An = ±1 transitions are forbidden, so we only plot the 
Rabi frequencies and detunings for the An = and An = ±2 transitions, and for a = 7r/4 the An = ±2 transitions 
are forbidden, so we only plot the Rabi frequencies and detunings for the An = and An = ±1 transitions. For these 
graphs the Lamb-Dicke parameter is taken to be 77 = 0.05, which is the value relevant for the experiments described 
in 0,[3ll,[32. 

V. RAMAN COUPLING FOR CESIUM 
A. Effective Hamiltonian 

So far we have discussed the FORT-Raman and Raman- Raman schemes in the context of a simple three-level model. 
In this section, we show how these schemes are modified when we take into account the multiplicity of levels in a 
physically realistic alkali atom, using cesium as an example. 

A level diagram for cesium is shown in Figure [H the levels relevant to our considerations include the ground state 
hyperfine manifolds 6S'i/2j^ = 3 and 6Si/2,F = 4, which correspond to ground states a and b of the three-level 
model, and the excited state manifolds 6P3/2 and 6P1/2, which correspond to the excited state e of the three- level 
model. The Hamiltonian for a free cesium atom is 

Ho = Y.''-\^^(^\ + \^HFiP4-Ps), (54) 

e 

where the sum is taken over all the states e in the QP3/2 and 6P1/2 excited state manifolds, and where P3 and P4 
are projection operators onto the F = 3 and F = A ground state manifolds. The quantity Ahf = (27r)(9.2 GHz) is 
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FIG. 3: Rabi frequencies and detunings for FORT wells with a = n/4: green curve is n„^„/f2p, blue curve is Qn^n+i/^p, 
red curve is 5„^„+i/(jj. 

the hyperfine splitting between the -F = 3 and F = A ground state manifolds, and LUe is the energy of excited state e, 
where the zero of energy is taken to be halfway between the two ground state manifolds. 

As in section fill Al we want to derive the Raman coupling that results when the atom is trapped within an optical 
cavity and one of the cavity modes is driven with a pair of beams that generate standing- wave fields inside the cavity. 
The coupling of the atom to the standing-wave fields is described by a Hamiltonian Hr that has the same form as 
equation (I14p . but with Q± given by 



^±-7(/±Aat)'/'V'(r)el-A 



(55) 



Here I± are the maximum intensities of the two fields, and 7 = (27r)(5.2MHz) and Isat = 2.19mW/cm^ are the 
spontaneous decay rate and saturation intensity for the 6P3/2 excited state manifold. We can express Igat £^s 



Isat = (4^V3)(7/AL), 



(56) 



where Xd2 — 852 nm is the wavelength of the 65*1/2 ~* 6P3/2 transition. The quantities e± are the polarizations of the 
two fields, which we will take to be linear and mutually orthogonal. Thus, the vectors {e+, e_ , k} form an orthonormal 
frame, where fc is a unit vector that lies along the cavity axis. The vector A is an atomic lowering operator, and is 
defined by 



^^^EEE E Pj'iF',F){F',m'\l,q;F,m)\6Pj,,F',m'){6S^/2,F,m\e;. 



(57) 
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FIG. 4; Level diagram for cesium. 



Here {F' ,m'\l,q;F,m) is the Clebsch-Gordan coefficient that connects ground state \F,m) to excited state \F',m') 
via polarization e* 



eo = z, e±i=T—^ix±iy), 



is a orthonormal basis of polarization vectors, and 



f3j>{F,F') = -(-1)^'V2J^TTV2^'TT| y^ ^/,^ "^' I , 



(58) 



(59) 



is a weighting factor for transitions between the 6Sit2,F and 6P,ji,F' hyperfine manifolds, where the quantity in 
brackets is a 6 J-symbol. 

Following the procedure we used in section IIII A[ we can adiabatically eliminate the excited states to obtain an 
effective Hamiltonian for the ground states. Because the derivation closely parallels the derivation given in section 
IIII Al we will omit the intermediate steps and simply quote the result: 



where 



He = -Ahf{P4 - P3) + Ve + {^e + n^) cos Snt, 



Ve - 5]^-^PF(f^+|e)(e|r2^ + A-|e)(e|f2U^F, 



F e 



^E = E^^3^+|e)(e|^-^4. 



(60) 

(61) 
(62) 



Here Ag = u^ — LOf. is the overall detuning of the Raman pair from excited state e. In the limit that Ag is much larger 
than the excited state hyperfine splittings, one can show that [33| 



^ ■^A,\e){e\A] = (1/3)[(2/Ad2 + 1/Aci)<5,y + 2i(l/Ai32 - 1/Ai5i)e„fe J^], 



(63) 



where A^i = wl — wdi and A£)2 = lol — tL)D2 are the overall detunings of the Raman pair from the cesium Dl and 
D2 lines. It is convenient to express these detunings as 



1/A 



Dl 



-(Al/2^)Cdi 



1/A 



D2 



\\lI2-k)Cd2. 



(64) 
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where 

Cdi = {\l/Xdi - l)-\ Cd2 = {Xl/\d2 - l)-^ (65) 

are dimensionless parameters. From equations (|6ip . (|62|) . and (|63p . we find that 

Ve = Ve IVX^)!', ^e = ^E IVXr)^ S, (66) 

where 

Ve = (7Vl2)((/++/-)//.at)(2/Ai32 + l/Am), (67) 

^E ^ {^V){I+I-/iI,)'/\1/Ad2 - l/Ani), (68) 



and 



E = 2P3fc-JP4 (69) 



is an atomic lowering operator that couples Zeeman states in F = 4 to Zeeman states in F = 3. If we collect 
these results, make the rotating wave approximation, and perform a unitary transformation to eliminate the time- 
dependence, we can express the effective Hamiltonian as 

He = -^(^4 - Ps) - Ve Hir^l' + ^^e |^(r-)p (S + ±^). (70) 

It is instructive to compare to the effective Hamiltonian for the three-level model given in equation (|24p with the 
effective Hamiltonian for the full cesium atom given in equation (j70p . The two Hamiltonians have similar forms, only 
the operator ax that coupled ground states a and b has been replaced by the operator S -|- E^ that couples Zeeman 
states within the ground state manifolds F — 3 and F = 4. In addition, we now have equations ((67)) and ([68)) that 
allow us to calculate the parameters Ve and Qe in terms of the intensities of the standing-wave fields. 

Following the same reasoning that was used in sections IIIICI and IIIIBi we can use the effective Hamiltonian given 
in equation ()70p to write down the total Hamiltonian H for the FORT-Raman and Raman-Raman configurations. In 
both cases the total Hamiltonian has the form H = Hext + Hint , where 

Hext = 7; hUp sin^ kFX, (71) 

Zm 



5 , , 1 

- 2(^4- Pa) +2^ 



H^nt = --{Pi- P:i) + -^p cos\kFX + a){t + t^). (72) 



fj2/2 n2/2 

Here Up — Upe^ '™^ and VLp — VLee^ '^ ''"^ are the axial trap depth and the effective Rabi frequency at radial 
coordinate p, and the parameters Up and U,e are calculated for the FORT-Raman and Raman-Raman configurations 
in the following sections IV Bl and IV CI 

B. FORT-Raman configuration 

As was discussed in section IlIIBI in the FORT-Raman configuration the FORT forms one leg of the Raman pair, 
and a weak Raman beam is added to form the second leg. The FORT resonantly drives mode np of the cavity, and 
the Raman beam drives the same mode as the FORT but is detuned from the cavity resonance by \5r\ ~ Ahf- 

We can obtain expressions for the FORT depth Up and the effective Rabi frequency ^e by using equations (|57)) 
and ([S5)) . which relate these quantities to the maximum intensities of the FORT and Raman beams inside the cavity, 
together with equation ()A5p from Appendix |^ which relates these maximum intensities to the optical powers of the 
FORT and Raman beams at the input of the cavity (note that because the Raman beam is detuned from the cavity 
resonance, its coupling into the cavity is suppressed). We find that 

Up = {-f/2AiT){^/Kp){2C^^ + C^^){Pp/P,), (73) 

^E - {ill2T:){^lnF){C^^-C^,){l + {2/^HF/nFfr''\PRPE/P^.Y'^. (74) 

Here Pp and Pr are the powers of the FORT and Raman beams at the input of the cavity, P^ is a reference power 
that is set by the cavity geometry and is defined in equation ()A4p of Appendix [^ Kp is the total energy decay rate for 
the FORT mode np, and C^2 ^nd G^^^ the detuning parameters at the FORT wavelength Af, are given by equation 
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([65)1 . It is interesting to note that for fixed powers in the FORT and Raman beams, the effective Rabi frequency i^E 
monotonically increases as the cavity decay rate np is reduced. 

In deriving the expression for Up given in equation ((73|) . we assumed that the detuning of the FORT from the 
cesium Dl and D2 hnes was the same for the F — 3 and F — 4 ground state hypcrfinc manifolds. This is a reasonable 
approximation, because these detunings are much larger than the hyperfine splitting Ahf- However, because the 
detuning of the F = 3 manifold is slightly larger than the detuning of the F — 4 manifold, the FORT potential is 
slightly weaker for F = 3. Thus, the FORT squeezes the two manifolds together, causing a small reduction in the 
effective hyperfine splitting. This effect, which is calculated in Appendix [Cl gives a slight position-dependence to the 
effective Raman detuning, but this can be neglected for many applications. 

C. Raman- Raman configuration 

As was discussed in section IIII C| in the Raman- Raman configuration the FORT resonantly drives mode np of the 
cavity, and pair of Raman beams drives mode nn of the cavity. We will assume that the two Raman beams have 
equal powers Pr and are tuned symmetrically about the cavity resonance. 

The FORT depth Up is given by equation (|73p . and we can obtain an expression for the Rabi frequency Qp by 
using equation (|68p . which relates the Rabi frequency to the maximum intensities of the Raman beams inside the 
cavity, together with equation (jASP from Appendix [XI which relates these maximum intensities to the optical powers 
of the Raman beams at the input of the cavity (note that because the Raman beams are detuned from the cavity 
resonance, their coupling into the cavity is suppressed). We find that 

np = (7/12vr)(7/«i^)(Ci^2 - CSi){l + {AHp/nj,f)-HPB./Pc)- (75) 

Here Pc is a reference power that is set by the cavity geometry and is defined in equation (jA4[) of Appendix El k^j 
is the total energy decay rate for the Raman mode nj^, and C|^2 ^^^ ^di' ^^^ detuning parameters at the Raman 
wavelength A/j, are given by equation (|65p . Note that for fixed powers in the Raman beams there is an optimal cavity 
decay rate kj^ = Ahp that maximizes the effective Rabi frequency flp. 

D. Zeeman transitions 

The operator E + S^ that appears in Hint couples individual Zeeman transitions between the F — 3 and F = 4 
ground state hyperfine manifolds. In this section, we calculate the matrix elements for these transitions. 

Let us introduce an arbitrary coordinate system {x, y, z} and define a set of Zeeman states {|3, m), |4, m)} relative 
to this coordinate system. We can express the unit vector k that lies along the cavity axis as 

k — cos (f) sin 9 x + sin sin 9 y + cos 9 z, (76) 

where 9 is the angle between the cavity axis k and the quantization axis z. Note that 

E = 2P3 fc • J P4 = P3 (2 Jz cos 9 + J+ e-"l' sin 9 + J- e"*' sin 9) F4, (77) 

where J± = JxHJy — =fV2^±i are angular momentum raising and lowering operators. Thus, the state |3, m) couples 
to states |4,to) and |4, ?n± 1), and the matrix elements corresponding to these transitions are 

{4,m\±''\3,m) = (1 - toV16)^/^ cos6', (78) 

{4,m+ 11^^3,771) = i(4 + m)i/2(5 + m)i/2e-^'^ sine, (79) 

8 

{4,m-l\t^3,m) = ^{4- fny/^{5 - m)^/^ e"^ sin9, (80) 

8 

where we have used that the matrix elements of J are given by [33| 

{F2,m2\J,\Fi,mi)=-{3/2)^/^{-lf'Vm + lly^ %l ^^H (^2, mall, g; Pi, mi). (81) 

Note that if the quantization axis is aligned along the cavity axis then Am = ±1 transitions are forbidden, and if the 
quantization axis is transverse to the cavity axis then Am = transitions are forbidden. 
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VI. RESOLVED-SIDEBAND COOLING 
A. Cooling schemes 

We have shown that the Raman couphng can drive transitions that raise or lower the axial vibrational quantum 
number n. In this section, we show how one can exploit these n-changing transitions to cool the axial motion to 
the vibrational ground state. We will first show how the cooling works using the three-level model, and then discuss 
cooling for a physically realistic cesium atom. 

One way to cool the atom is to alternate coherent Raman pulses tuned to an n-lowering transition with incoherent 
repumping pulses. To see how this works, let us assume that we start out with the atom in state \a, n). We can lower 
the vibrational quantum number by driving the atom with a coherent Raman pulse tuned to the n — > n— 1 transition, 
which transfers some of the population from |a, n) to |6, n — 1). The atom can then be repumped to ground state a 
by driving the h — e transition with near-resonant light. The repumping light drives the atom to the excited state, 
from which it spontaneously decays to either ground state 6, where it continues to interact with the repumping light, 
or to ground state a, where it is dark to the light. If the atom is sufficiently cold to begin with, then the repumping 
process is unlikely to change the atom's vibrational state, because the matrix elements for n-changing transitions are 
suppressed relative to the matrix elements for n-preserving transitions by at least r]f,y/n, where rje = {2muj)~^''^ uje- 
Thus, the net effect of the Raman and repumping pulses is to move some of the population from state |a, n) to state 
\a,n — 1). By iterating the pulse sequence, the atom can be cooled to a state that has a mean vibrational quantum 
number n that is close to zero. 

The same type of scheme can be used to cool a multi-level alkali atom. For a cesium atom, the F — 3 ground state 
manifold plays the role of state a and the F = 4 ground state manifold plays the role of state b: we start with the 
atom in a random Zeeman state in i^ = 3, drive the atom with a coherent Raman pulse tuned to the n —f n — 1 
transition, and then repump the atom to _F = 3. It is easiest to understand the effects of these pulses if we choose the 
quantization axis to lie along the cavity axis, so that only Am = transitions are allowed and the Raman coupling 
drives transitions between pairs of states |3, m) ^^ |4, m). If the ambient magnetic fields are nulled, then these Zeeman 
transitions are all degenerate, so the transition frequency of the n —f n — 1 transition is independent of m. Thus, 
the coherent Raman pulse is effective at lowering the vibrational quantum number regardless of which Zeeman state 
in _F = 3 the atom started in: each Zeeman pair behaves equivalently, except for a slight m-dependence in the Rabi 
frequency that is given by equation (|78p . The Zeeman state of the atom is then scrambled during the repumping 
phase, so at the beginning of the next cooling cycle the atom starts out in a potentially new Zeeman state. 

The amount of time it takes to cool the atom is determined by the amount of time it takes to repump the atom, 
which is set by the spontaneous decay rate of the excited state, and by the amount of time it takes to perform the 
coherent Raman pulse, which is set by the Rabi frequency i^n->n-i- The cooling rate can be increased by increasing 
the Rabi frequency, but as we increase the Rabi frequency we begin to off-resonantly drive the n —^ n transition, and 
this sets an upper limit to the Rabi frequency that can be used. Off-resonant driving of the n — > n transition becomes 
important when Qn-^n ^ ^, so the upper limit to the Rabi frequency is given by f2„^„_i ^ rj^/nuj. 

There is also a lower limit to the value of n that can be achieved with this cooling scheme, which is set by two 
different factors. First, when we resonantly drive the n ^ n — 1 transition with the coherent Raman pulse, we can 
also off-resonantly drive the n — > n -I- 1 transition. This mechanism gives a lower limit oi n ^ (rio-+i/2(j-')^. We 
can reduce this limit by reducing the Rabi frequency, but since the Rabi frequency determines the cooling rate, this 
also slows down the cooling. In addition, there are problems with using small Rabi frequencies that are due to the 
anharmonicity of the FORT, which will be discussed later. Ideally, one would gradually reduce the Rabi frequency as 
the atom cools, so as to balance the conflicting demands for a high cooling rate and a low value of n. A second factor 
that limits n is the fact that when the atom is repumped it will not always remain in the same vibrational state it 
started out in, since the Lamb-Dicke suppression of n-changing transitions is not perfect. This mechanism gives a 
lower limit of n ~ 77^ . 

The cooling scheme described above can be modified in several ways. First, rather than alternating Raman pulses 
with repumping pulses, it is also possible to continuously drive the atom with both Raman and repumping light, 
and this is the method that was used in 30] and [31]. Second, the cooling scheme we described relies on An = —1 
transitions, but it is also possible to cool the atom using An = —2 transitions. Indeed, for the FORT-Raman 
configuration the An = —1 transitions are forbidden, so the atom can only be cooled using An = —2 transitions. 
Cooling via An = —2 transitions tends to be slower than cooling via An = — 1 transitions, since the condition 
f2„_>„ ~ uj gives an upper limit on the Rabi frequency of Q,n—>n-2 ^ rfnuj. Also, for An — —2 transitions both the 
state ]a, 0) and the state ]a, 1) decouple from the Raman pulse, so the state to which the system cools depends on the 
initial state: if we start in a state ]a, n) with n even then the system cools to ]a, 0), and if we start in state ]a, n) with 
n odd then the system cools to ]a, 1). 

Note that because the FORT is anharmonic, the resonant frequency of the An = — 1 and An = —2 transitions 
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depends on the value of n. This means that if we keep the Raman detuning S set at at a fixed value throughout 
the cooling process, then the detuning of the Raman pulse from the atom will change as the atom cools. We can 
estimate the importance of this effect by considering An = —1 and An = —2 transitions as separate cases. First we 
will consider An = — 1 transitions. Let us assume that we set the Raman detuning to S = —uj, so the detuning of the 
Raman pulse from the n — > n — 1 transition is 

A„^„_i =5 - 5n^n-i ^ rfnuj. (82) 

As we have shown, the maximum Rabi frequency that can be used is J7„^„_i ^ rj^/nuj^ and for this maximum 
value the ratio of the detuning to the Rabi frequency is ~ Tj^/n^ which is small for cold atoms. Thus, for An = — 1 
transitions the dependence of the detuning on n is a small effect; we can simply set the Raman detuning to 5 — — w, 
and as long as the atoms start out reasonably cold the cooling will always be efficient. 

Now consider An = — 2 transitions. We will assume that the Raman detuning is set to (5 = — 2w, so the detuning 
of the Raman pulse from the n ^ n — 2 transition is 

A„^„_2 ^ 5 - &n^n~2 ~ ^rj^nuj. (83) 

As we have shown, the maximum Rabi frequency that can be used is J7„^„_2 ^ rj'^nuj, and for this maximum value 
the ratio of the detuning to the Rabi frequency is ~ 2. Thus, for An = —2 transitions the dependence of the detuning 
on n is a significant effect. To compensate for this problem, one could slowly decrease the Raman detuning 6 during 
the cooling process to ensure that the Raman pulse remains in resonance as the atom cools. 

Although we have focused on cooling the axial motion of the atom, it is possible to implement the axial cooling 
schemes in such a way that they cool the atom's radial motion as well. This is accomplished by configuring the 
repumping light so that it provides polarization gradient cooling [34] in the plane transverse to the cavity axis. 
Specifically, the repumping light is tuned blue of the 6S1/2, F — 4 ^ 6P3/2, -F" — 4 transition, and is delivered to the 
atom via two pairs of counter-propagating beams. The two pairs of beams are perpendicular to one another and to 
the cavity axis, and therefore provide cooling in both transverse directions. 

B. Measuring the temperature 

One can characterize the effectiveness of the cooling schemes described in the previous section by using Raman 
spectroscopy to measure the temperature of the atom. In what follows, we will assume that An = — 1 transitions are 
used to cool the atom, although the same methods can also be applied to cooling via An = — 2 transitions. 

To measure a Raman spectrum, we cool the atom, pump it into ground state a, and then drive it with a coherent 
Raman pulse with detuning 6. We then check if the atom was transfered to b. By iterating this sequence one can 
measure the probability that the atom is transfered from a to 6 by the Raman pulse, and by repeating this measurement 
for Raman pulses of different detunings one can map out a Raman spectrum. For an atom in vibrational state n, the 
Raman spectrum will exhibit a peak at (5 = 0, which corresponds to n — > n transitions, and peaks at i5 = ±(5„, which 
correspond to n — > n ± 1 transitions. We will refer to the peak at (5 = as the carrier, and the peaks at i5 = — 15„ and 
(5 = (5„ as the red and blue sidebands. Because of the FORT anharmonicity (5„ depends on n, but we will assume that 
the atoms are cold enough that this effect can be neglected and simply take (5„ ~ lj. 

One way to determine the axial temperature of the atom is to measure the ratio of the red to the blue sideband; this 
is the same technique as was used in [25| to determine the temperature of a trapped ion. For a thermal distribution, 
the probability that the atom has axial vibrational quantum number n is given by 

where n = (e^'^ — 1)~^ is the mean vibrational quantum number. If we start with the atom in state a and resonantly 
drive the blue sideband with a Raman pulse of duration t, the probability that the atom is transfered to state b is 
given by 



Pb 



^P„sin2(r!„_„+it/2). (85) 



If we start in state a and resonantly drive the red sideband, the probability that the atom is transfered to state b is 
given by 

00 00 

Pr = J2^^ sin^{^n^n-it/2) = ^ P„+i sin^ (a„+i^„i/2) . (86) 

n=l n—0 
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Note that 

so the ratio of the transfer probabihties for the red and blue sidebands is 

Pr/Pb^nl{n + l). (88) 

An alternative way to quantify the cooling is to measure the population in the vibrational ground state. This can 
be accomplished by pumping the atom to state a and applying a Raman pulse whose detuning is adiabatically swept 
across the red sideband. If the atom started in the vibrational ground state n — Q then it will remain in state a, and 
if the atom started in a vibrational state n > then the Raman pulse will adiabatically transfer it to state h. Thus, 
the population in the vibrational ground state is given by the probability that the atom remains in state a after the 
adiabatic sweep has been completed. The advantage of this method is that it does not rely on assuming that the 
atoms are thermally distributed. 

It is also possible to use Raman spectroscopy to say something about the radial temperature: since the axial 
frequency Lo{p) = LOa e"'' '""^ depends on the radial coordinate p, the width of the sidebands depends on the radial 
temperature. The probability that the atom has axial frequency A is given by 

1 P°° 
P{^)^^j e-P^^PU{A-u;{p))pdp, (89) 

where U{p) ~ —Upe''^^ /^^ is the potential for radial motion, and 

/•OO 

Z= / e-^^'^P^ pdp. (90) 

Jo 

If the radial temperature is small compared to the trap depth {13Uf ^1), then we can make a harmonic approximation 
and perform the integral analytically: 

p(A) = {l/uja){2pUF) e{\ - A/wa) exp(-2/3C/j.(l - A/c^,)). (91) 

Thus, if the blue sideband has width 5lj^ one can put an upper limit on the radial temperature of fcsT < 2Uf{5uj /lOo). 

C. Cooling simulation 

The cooling schemes discussed in section IVI Al can be simulated on a computer. We will take the Hamiltonian for 
the system to be 

H = Hint + Hext + Hl, (92) 

where Hint and H^^t ai'e given by equations (j39p and psp . and where 

Hl = - Ap |e>(e| + (f^p/2) (|6) (e| + |e) (6|) (93) 

describes the coupling of the atom to repumping light. Here fip is the Rabi frequency of the repumping light and 
Ap is the detuning of the light from the h — e transition. As was discussed in section IVI Al in order to radially cool 
the atom we use light that is blue-detuned from the 4 — 4' transition to repump the atom. To model this in the 
simulation, we will assume that the excited state decays to ground state a at rate "fa = (5/12)7 and to ground state 
h at rate 7;, = (7/12)7, where 7 — (27r)(5.2MHz) is the spontaneous decay rate for the QPz/2 manifold of cesium, and 
the prefactors 5/12 and 7/12, the branching ratios for spontaneous decay on the QP^/2tP = 4 ^ 6S'i/2,^ = 3 and 
6-F3/2J F = A ^r 6S'i/2, F — A transitions, are given by equation ((59| . Also, we will take Ap = (27r)(10MIIz) to be the 
detuning that optimizes the polarization-gradient cooling. 

We can write down a master equation for the system, which describes both the coherent evolution due to H and the 
incoherent evolution due to the spontaneous decay from the excited state. Given an initial state, we can numerically 
integrate the master equation to obtain the state of the system at later times. In Figure [5^, we use this method 
to simulate cooling in the Raman- Raman configuration: we start the system in state |a, 5) and plot n as a function 
of time. We assume that the atom is trapped in a FORT well with a — it /A, and use An = — 1 transitions to 
perform the cooling. For this simulation, the cooling parameters are Vlp — (27r)(0.2MIIz), Vlp — — (27r)(3MHz), 
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FIG. 5: Cooling simulation: mean vibrational quantum number n versus time t, (a) cooling via An = — 1 transitions starting 
from state |a,5), (b) cooling via An — —2 transitions starting from states |a, 3) and |a, 4). 
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FIG. 6: Cooling simulation: (a) asymptotic n versus fip, (b) asymptotic n versus Qp, (c) asymptotic n versus 5. 



6 = — (27r)(0.5MHz). In Figure [Sb, we simulate cooling in the FORT-Raman configuration. In the FORT-Raman 
configuration a = for all the FORT wells and the An = — 1 transitions are forbidden, so we use An = —2 transitions 
to perform the cooling. As was previously discussed, this means that the asymptotic state to which the system cools 
depends on the initial state. Two curves are shown in the graph: for one, we start the system in state |a, 3); for 
the other, we start the system in state |a,4). For these simulations, the cooling parameters are fip — (27r)(0.2MIIz), 
rip = -(27r)(3MHz), S = -(27r)(0.75MHz). 

In addition to simulating the time-evolution of the system, we can calculate the asymptotic value of n by solving 
the master equation for the steady-state density matrix. This can be used to study the dependence of the asymptotic 
value of n on the various cooling parameters. In Figure [S] we consider cooling in the Raman- Raman configuration 
for atoms with 7r/4, and plot the asymptotic value of n as a function of f2p, Op, and d. The parameters that are not 
being varied are set to the same values used for the cooling simulation shown in Figure [5^. These graphs show that 
the cooling scheme is quite robust and works efficiently over a broad range of parameters. 



VII. CONCLUSION 



We have described two schemes for driving Raman transitions in an atom trapped within a high-finesse optical 
cavity. These schemes can be used to control both the internal and motional degrees of freedom of the atom, and 
provide powerful tools for studying cavity QED; as an example, we have shown in detail how the Raman schemes 
can be used to cool the atom to the quantum ground state of the trapping potential. Although the two schemes are 
similar in many respects, they do have some important differences. The FORT-Raman scheme has the advantage 
that the Raman coupling is independent of the FORT well in which the atom is trapped, and is thus better suited 
for manipulating the internal state of the atom. On the other hand, the Raman-Raman scheme has the advantage 
that the n ^ n ± 1 transitions are allowed for most FORT wells, and is thus better suited for cooling. The ability to 
coherently control the atom is a key requirement for many cavity QED protocols, and these Raman schemes should 
open up new possibilities for experiments in cavity QED. 
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APPENDIX A: CAVITY MODE STRUCTURE 

Here we describe the mode structure of the optical cavity. The cavity we will be considering consists of two 
symmetric mirrors of radius R that are separated by a distance L. It is convenient to define a cylindrical coordinate 
system that is centered on the cavity axis: we will denote the distance from the cavity axis by p and the displacement 
along the cavity axis by z, where the mirrors are located at z = and z — L. The cavity supports a set of discrete 
modes with resonant frequencies at integer multiples of the free spectral range vfsr = l/2i, where for each frequency 
there are two degenerate modes corresponding to the two polarization states transverse to the cavity axis. Consider 
one of the polarization modes with mode order n. We will let uj = 2Tmi>FSR. denote the resonant frequency of the 
mode, and let A = 2'k/uj and k = 2tt/X denote the corresponding wavelength and wavenumber. We can characterize 
the shape of the mode by a dimensionless function ip{'r) that is given by 

tlj{f}=e-P^^'"\mkz, (Al) 

where w — {L{2R — L)/k'^)^'* is the mode radius. If we drive the cavity with an input beam that has power P^ and 
frequency Ui, then the intensity at a point f inside the cavity is 

/(f) =. (2/a.T/)(1 + {2A/Kf)-' mr)\^ P^, (A2) 

where A = uji — lUc is the detuning of the input beam from the cavity resonance, k is the total energy decay rate for 
the mode, and V^ the mode volume, is given by 

V= f |V'(f)p dV = (Ai/8)(2i?L)i/2. (A3) 

In order to relate the input power to the maximum intensity inside the cavity, it is convenient to define a power 

P, = {2V/\)hat = (i/4) (2i?i) 1/2 /^^^ . ( A4) 

Note that Pc is the same for all cavity modes; it depends only on the cavity geometry, not the mode number. We can 
then express the maximum intensity inside the cavity as 

Ima. = (kA)-1 (1 + (2A/At)2)-1 (P,/P,) /,,,. (A5) 

APPENDIX B: DISTORTION OF THE TRAPPING POTENTIAL 

In the Raman-Raman configuration, the lack of registration between the FORT and Raman beams causes a well- 
dependent distortion of the trapping potential. Here we calculate this effect. From equation ([32)l. we see that the 
total potential for the Raman-Raman configuration is given by 

U{r) = -Uf e-^P^/""^" sin^ kFZ ~ Ur c-^pV"-?, gin^ fc^^. (Bl) 

We will assume that an atom is trapped in well r of the FORT, and define a coordinate x = z — Zr and a phase 
c^ = (fc-R ~ kF)zr. Note that 

sin^ kFZ — cos^ kFX, siri^ kfiz — cos'^{kiix + a). (B2) 

We will assume that the FORT and Raman beams drive nearby modes of the cavity, so {kn — kF\ <^ k^, kF ■ In this 
limit, we can approximate U{r) by replacing wr with wf and replacing k^x with kFx: 

U{r) = ~Uf e-'^P^I-^^" cos^ Ufx ~ Ur g-^pV^-J cos^{kFX + a). (B3) 
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As in section ini we will take the atom to be radially stationary, and take p as a constant parameter that enters into 
the potential for axial motion. We can then write the potential as 

U{r) = Up sin^ikFX + 9), (B4) 

r) 2 / 2 

where Up = Uq e~ '' '""^ is the axial trap depth at radial coordinate p. The quantity Uq is given by 

Uo = {U^ + 2UFURCOs2a + U]i)^^\ (B5) 

and 9 is given by 

tan 20 = {Uf + L/i^cos2a)~^[/fl,sin2a. (B6) 

APPENDIX C: DIFFERENTIAL STARK SHIFT 

The FORT potential is slightly weaker for the F = 3 ground state manifold than for the F — 4 ground state 
manifold, so there is a small differential Stark shift. Here we calculate this effect. The Hamiltonian that describes 
the differential Stark shift is 

Hd = ^Sd \^Fir)\'' {Pi - Ps) =Hd = ^6p (P4 - P3) cos^ kFX, (CI) 

where 5d is the differential Stark shift at an intensity maximum and 5p = 5d e~^^ ' "'^ is the maximum differential 
Stark shift at radial position p. We can calculate 5f, as follows. The FORT depth Uf is given by (p7|) with /+ = If 
and /_ = 0: 

Uf = (7Vl2)(/F//.at)(2/A|;2 + l/^^,). (C2) 

Thus, the differential Stark shift at an intensity maximum is given by 

j^ If f 2 1 \ j^ If f 2 1 



Expanding to first order in A//^, we find that 

dn = -Uf i2{C^o2? + (C^i)')(2C^2 + C^diY' (Ahf/u^f)- (C4) 

where C^i and C|^2i the detuning parameters at the FORT wavelength Xf, are given by equation (|65p . 
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